Anti-localization of graphene under substrate electric field 
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A simple criterion is provided iiow tiie (anti-)localization properties of grapliene are determined 
in the presence of inter-valley scattering, Kane-Mele topological mass term, and Rashba spin-orbit 
interaction (SOI). A set of (pseudo) time-reversal operations show that the number of effective 
internal degrees of freedom, such as spin and pseudo-spins distinguishing the sublattice and the 
valley, is the crucial parameter for localization. It is predicted that perpendicular electric field due 
to gate voltage of the substrate drives the system to anti-localization by enhancing the Rashba SOI. 



Graphene has a strong tendency not to localize 
reflecting its linear spectrum [j]. In ordinary two- 
dimensional metals, scaling either to weak localization 
(WL) or to weak anti-localization (AL) is controlled by 
impurity scattering with spin-orbit interaction If 
time reversal symmetry (TRS) is broken by external or 
internal magnetic field, the system exhibits neither WL 
nor AL , and belongs to the unitary class 0] . Graphene 
shows the three localization classes even without mag- 
netic impurities, but with potential scatterers. 

In the literature, theory predicts that graphene under 
doping shows AL provided inter-valley scattering can be 
neglected [H, 0, S| . Together with the absence of back- 
ward scattering, AL in graphene is a clear manifestation 
of Berry phase tt 0] . Inter- valley scattering, on the other 
hand, drives the system from AL to WL Experi- 
ments [1, show also a unitary behavior. Absence of 
WL may be attributed to ripples [ll|, It is natu- 

ral to ask how the localization properties are influenced 
by modification of the linear spectrum by a finite mass, 
whose magnitude is under debate in pho toemission ex- 
periments [H, and theoretically [Tsl . ITg} . Another 
modification is due to the Rashba SOI which inevitably 
appears in graphencs on substrates. Kane and Mele have 
proposed that graphene with intrinsic spin-orbit interac- 
tion realizes a Z2 topological insulator [l^. This paper 
demonstrates another striking aspect of graphene under 
doping that crossovers between different symmetry class 




FIG. 1: Weak localization classes of graphene-based mod- 
els with different types of the mass and impurity scattering. 
WL, AL and U refer to weak localization (orthogonal class), 
weak anti-localization (symplectic class) and absence of WL 
(unitary class), respectively. 



of localization are controlled simply by the number of 
effective internal degrees of freedom. 

The inter-valley coupling depends on the range of im- 
purity potential 17]. Long-range scatterers (LRS) do 
not involve intervalley scattering, whereas short-range 
scatterers (SRS) activates the valley spin. Localization 
properties of a disordered Kane-Mele model have been 
studied numerically [H, [l^ . Since tight-binding approx- 
imation is employed there, inter-valley scattering is al- 
ways present, and cannot be controlled. The strategy of 
our study is from the opposite direction; we employ the 
lowest-order WL theory, but analyze systematically each 
element that influences the localization behavior. 

Our results are summarized in FIG.l, which shows lo- 
calization classes WL, AL, and U depending on the inter- 
valley scattering, types of the mass term, and the Rashba 
SOI. Without these effects, the system becomes massless 
and belongs to AL, as is the case of single Dirac cone 
[l|,[l]. Note that, in the presence of inter- valley scatter- 
ing, crossovers to AL occurs by switching on the Rashba 
SOI. The crossovers should be achieved by controlling 
the gate voltage of the substrate on which the system is 
placed. In the single valley case, the Rashba SOI leads to 
either WL or U, but not to AL. Such single valley system 
appears in HgTe/CdTe heterostructure |20l. [2lj|. 

We take the doped and disordered Kane-Mele model 
for our weak-coupling perturbative theory. The Kane- 
Mele model has a three story structure, depicted in TA- 
BLE I: (i) graphene in the massless limit, (ii) topolog- 
ical mass term, encoding Kane-Mele spin-orbit interac- 
tion, (iii) Rashba term, playing the role of activating the 
real spin degree of freedom. Note that the Kane-Mele 
model has also valley degree of freedom, corresponding 
to two Dirac points of graphene: K and K' . The ionic 
mass term is induced by chiral symmetry breaking stag- 
gered chemical potential. The Kane-Mele model pos- 
sesses three types of pseudo or real spins, represented 
by Pauli's matrices, (?, r and s, operating in different 
subspaces: a acts on the sublattice spin A-B, r on the 
valley spin K-K' , and s on the real spin. In the contin- 
uum limit, the Kane-Mele Hamiltonian, 

Hkm = Hi+Ha + Hr, (1) 

consists of the following three elements: (i) Hi = 
hvpipxCTxTz +PyO-y), (ii) Ha = -Aa^r^^Sz, (iii) Hr = 
— ^nicxTzSy — aySx)/2, each describing the correspond- 
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TABLE I: Three story structure of Kane-Mele Z2 topological insulator, and its WL properties under doping. LRS is equivalent 
to the single valley model. The parity of Ng^, the number of activated spin degrees of freedom, determines its WL properties: 
standard WL (orthogonal) or AL (symplectic). Broken TRS leads to U (unitary) behavior. 





LRS (single valley) 


SRS {K-K' coupled) 


(i) massless graphene: }i\ = PxCrrTz + PyUy 


Ns = 1 {AB) AL 


TV, = 2 [AB,KK') WL 


(ii) mass terms; H2 = -ffi + HA,m 

(a) topological —AazT^s^ vs. (b) ionic ma^ 


unitary 

no l/g-correction 


(a) unitary 

(b) Ns^2 {AB,KK') WL 


(iii) Rashba spin-orbit interaction: 

H3 = H2 — ^RiOxTzSy — (JySx)/2 


(a) Ns = 2 {AB, real spin) WL 

(b) unitary 


iV, = 3 [AB, KK', real spin) ^ AL 



ing floor of the three story structure. For comparison 
with i?A, we consider also ionic mass term Hm = ma^. 
In the ionic mass case, contributions to a^y from the K- 
and if'- valleys cancel [l^l • In the topological mass case, 
this cancellation of anomaly does not occur any longer 
[2^ , and a quantized spin Hall effect with preserved TRS 
occurs. In the absence of Rashba term Hr, the Hamilto- 
nian Hkm is diagonal in the real spin s space, implying 
that the latter is actually inactive. 

As shown in TABLE I, sublattice spin a is always ac- 
tive, whereas valley and real spins can be switched on 
and off, leading to four different cases for the subspace 
S spanned by active spins. In order to distinguish active 
and inactive spins, we introduce (pseudo) TRS opera- 
tions Ts, defined in the subspace E, such that 



(2) 



where S — {<?}, {(?, r}, {a,s}, {a,T,s}. Their explicit 
forms are given by 



Tfj — i<7yC, Tar — ^xC , 



Tx{-'iSy)C, 



where C is complex conjugation. T„rs represents the 
genuine TRS operation. Effective TRS of the system is, 
therefore, determined by the transformation property of 
the mass term (see TABLE II). When a mass term is 
odd against TRS, the system shows the unitary behav- 
ior. Four unitary phases in FIG.l correspond to the four 
minus signs in TABLE II. If some (pseudo or genuine) 
TRS exists in the system, its weak localization prop- 
erty is determined by the number Ng of the activated 
spin degrees of freedom. One can verify = 1 if Ng is 
even, whereas T| = — 1 if A^s is odd. The former (latter) 
corresponds to the orthogonal (symplectic) class in the 
random matrix theory Q , and leads to constructive (de- 
structive) interference between two scattering processes 
transformed from one to the other by T^. 

In the presence of a mass term, irrespective of its type, 
one finds unitary behavior for LRS. Scattering matrix 
elements are diagonal in r-space, and a is the only active 
spin {Ns = 1, TRS broken). In the massless case, a 
pseudo TRS operation T„ mimics the role of genuine TRS 
|24| . Once TRS is effectively restored, the system's WL 



property is determined by the parity of Ns- The mass 
term, on the other hand, explicitly breaks T„. 

SRS activate the valley spin f {Ns — 2), since its ma- 
trix elements involve off-diagonal terms in this subspace. 
In the case of ionic mass term, this leads the system to 
standard WL, since activation of the valley spin f re- 
stores the (pseudo) TRS. In the case of topological mass 
term, the system stays unitary, since restoration of TRS 
needs also the activation of real spin s. The latter is em- 
bodied by the Rashba SOI. In the presence of both SRS 
and Rashba SOI, we predict AL, since Ng = 3. 

Starting with the case of Aj,; = 0, let us go into some 
details of diagrammatic calculations. In the presence of 
topological mass term H/^, The spinor part of the con- 
duction band eigenstates reads. 



\Ka) 
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, \K'a) = 


e"^° sin ^ 






gi0„ sin ^ 




— cos ^ 



(3) 



where a specifies a three-dimensional fictitious momen- 
tum p = (p,r,p.y, — A), and the polar angles 9, (j) satisfy 

cosO = — A/^ 



Px+Py 



COS( 



■Px/^x 



■pI- 



LRS have a potential range much larger than the inter- 
atomic distance, and do not couple K and K' . The sys- 
tem cannot see the difference between two types of mass 
term, both showing unitary behavior, i.e., the diffusion 
type singularity is cut off by a cooperon's lifetime. This 
unitary phase shows a crossover to the well-established 
syrnplectic behavior of graphene in the single Dirac cone 



TABLE II: Time reversal operations Te, relevant in the sub- 
space spanned by activated spins. Transformation property of 
a mass term O = maz, l^azTzSz under Ts: T^OT^^ = ±0. 
The sign appears in the table. U refers to unitary class. 
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FIG. 2: Particle-particle ladders, (a) Relevant diagrams 
in the presence of short-range scatterers (SRS). "cis" and 
"trans" refers to specific configurations of the valleys, (b-d) 
Bare diagrams involving inter-branch processes {Xr 7^ 0, no 
SRS). (b) 7_+ contributes to 1/g^-singularity, whereas such 
diagrams as (c) and (d) are irrelevant to the singularity, since 
Pa + Pa' cannot be smaller than the order of Xr. 



SRS involve inter- valley scattering, allowing for distin- 
guishing the two different types of mass term: topological 
and ionic. The scattering matrix elements involve a pro- 
jection operator in the AB sublattice space, 'Pa,b- As for 
the singular contribution, one can focus on the two types 
of diagrams shown in FIG.l. The "trans" component 7* 
reads explicitly. 



7t = 2Tn^nAu\{K' PIVat 
+ 2^:vnBu\{K'|3\VB^ 



= -e'('^°"'^'')77ssin2 6l/2 (= 



\Ka){Kp'\VAT+\K'a') 
\Ka){K(3'\rBT+\K'a') 

7c), (4) 



')/2, 



where we have introduced rjs — 2Tiv{nAUA + nBUg 
with V, nA,B and ua.b being, respectively, the density 
of states, the impurity density and the typical strength 
of scattering potential at the A (B) sites. 74 has an ad- 
ditional minus sign, which plays the role of driving the 
crossover from symplectic to orthogonal behavior in the 
massless limit The Bethe-Salpeter equation (BSE) 
takes the form of two coupled equations: 



J Q/3 



7c 
7t 



a/3 



7c It 
It Ic 



(5) 



After diagonalization, one finds, Tc 
l-(7 



(1) 



h Tf = 0, and 
7i, where Us — 



^ (r, - Vt) = 7. 

rs(l — TsDq^) with T5 — 1/r/s being the scattering time. 
We also introduced 7^^' in the light of a general expres- 
sion: 7 = 7(')g«'(0Q-00)_ Cancellation between T5 and 
the bare 7's is incomplete, giving a finite lifetime. The 
system shows a unitary behavior, irrespective of the pre- 
served TRS of underlying Hamiltonian H = Hi + Ha- 
1 /q^-singularity is recovered in the limit E ^ 00 (orthog- 
onal class). 

In the case of ionic mass term, a complete cancellation 
between the self-energy and the bare vertex function oc- 



curs, driving a crossover from unitary to orthogonal sym- 
metry class. This situation is relevant to the intrinsic sin- 
gle valley system [l^, HH in the presence of off-diagonal 
interaction. 

Rashba SOI appears when inversion symmetry with 
respect to the 2D plane is broken, say, by a perpen- 
dicular electric field. Rashba SOI lifts the two-fold 
real spin degeneracy of the two conduction and two va- 
lence bands. An accidental degeneracy occurs on top of 
the valence bands, whereas the conduction bands Eu± 

are spht by 2Ai^: E± = 



Px 



p2 + (A±Afl/2)2 ± 



\r/2. In order to parametrize the corresponding eigen 
spinors, it is convenient to introduce fictitious 3D mo- 
menta: p± = {px,Py, ^ ± Xr/2), and cos9± = (A ± 



Afl/2)/y/^ + p2 + (A ± Afl,/2)2. In terms of these pa- 
rameters, the eigen spinors corresponding to E± read. 



\K±) 



cos 



sm 
1 



=F« cos ^ 
Tie"l' sin ^ 



\K'±) 



cos -f- 
-e-">' sin% 
Tie"^ sin ^ 
±1 cos 



(6) 

where we omitted the normalization factor 1 / V2. Rashba 
coupling imposes a stronger constraint on the choice of 
our basis. As a result, the matrix element, such as, 
(K(3\l\Ka) = cos{(t)a ~ sin^{6^/2) + cos2(0_/2), be- 
comes real (no Berry phase). 

When E > A + A_r, inter-branch matrix elements plays 
a role. They modify the scattering time t± for the \Kuzt) 
branches as. 



1 



2 



sm h 2 cos — 



sm 



1± 
2 



(7) 



As for particle-particle ladders, such as FIG. 2 (b), four 
electron states a, f3, a', f3' can, in principle, take either of 
the two channel indices, \K—) or \K+), generating eight 
types of diagrams in total. However, a simplification is 
possible at this level: since we are interested only in the 
l/g'^-singular part of cooperon diagrams, we need Pa + 
Pa' — q 0. This means that a and a' must belong 
to the same branch. We can thus safely focus on such 

diagrams as 7-1 and 7 |- depicted in FIG. 2 (b), or even 

simpler 7__ and 7++. Other 7's such as FIG. 2 (c) and 
(d) are irrelevant to singularity. Explicit form of 

relevant 7's are given as, 



7±± = VL 



cos 



2cos^ 



7±T = -VLSm 



, - 0/3) sin^ — 
'^/3)sin —cos — - 
1 - (/),3)sm — sm— . 



cos 



0± 



(8) 



These four types of diagrams satisfy coupled BSE. But 
its 4 X 4 coupling matrix is shown to be block diagonal- 
izable. One can, e.g., decouple F and r_| , from the 
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remaining part. The former two obey coupled equations, 
which take symboHcally the following form: 
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1-- 


n_ 7_+n+ 




















n- 7++n+_ 






(9) 


where Il-t 




T±(l - 




with T± given 


in Eq.©. 



To identify the singular contribution, we employ the ex- 
pansion into different angular momentum contributions, 
both for 7's and T's as T = r(')e*'('^°-'^''\ and pick 
up only the / = component. The following identity, 



det 



1 



(0) 

77_T_ 



(0) 



(0) 1 (0) 



(10) 



ensures that the dressed cooperons show indeed l/q^- 
singularity at the I = channel. Rashba SOI thus drives 
the system to standard WL, whenever the Fermi level is 
above the gap. 

If one adiabatically switches off the Rashba term, the 
simplification we have made for justifying Eq.® is no 
longer valid. In the limit of vanishing Ai? we cannot sim- 
ply neglect such diagrams as FIG.l (c), (d). They could 
contribute equally to the 1 /q^ singularity if the singular- 
ity ever appears. Relations such as Pa + Pa' — q — can 
be satisfied in these diagrams. In the coupled BSE, the 
cooperons acquire more channels to couple to, and one 
can no longer decouple F and F_| as Eq.®. As a re- 
sult, the cancellation property between particle-particle 
ladders and the self-energy such as Eq. fTUI) is lost. The 
loss of cancellation property immediately leads to the ab- 
sence of WL. On the contrary, if one starts from the situ- 



ation where Rashba SOI is absent in the first place, it is 
possible to choose a basis in which that two (degenerate) 
channels are decoupled. 

As for the ionic mass, Rashba SOI plays no role and 
the system stays unitary. The same applies to the case 
of ripples J^], the latter possessing the same symmetry 
properties as the ionic mass. 

We finally consider the most realistic case of SRS in 
the presence of Rashba SOI ( A_r 7^ 0). One has to use 
([5]) instead of ([3]), and develop a similar analysis. One 
finds in this case the cancellation between the self-energy 
and 7 at the I = channel. But here additional mi- 
nus sign analogous to Eq.Q appears due to inter- valley 
scattering, driving the system to AL. 

Let us estimate the strength of gate electric field 
required for observing the crossover to AL. For the 
crossover to be experimentally accessible, Rashba SOI 
needs to be the order of ~ 1 K. This corresponds to the 
electric field of order ~ IV/nm il5il . a value attainable in 
double-gated graphene devices [2a| • The crossover to AL 
will be observed for a sample with insignificant ripples. 
A similar crossover due to Rashba SOI has been observed 
in another context in InGaAs/InAlAs quantum well [26j . 
In the case of graphene, crossovers occur not only from 
WL, but also from the unitary class [§]. 

In conclusion, we have identified key elements respon- 
sible for rich localization properties of graphene that are 
controllable by substrate electric field. The criterion 
whether the system is driven to WL or AL is simply 
given by the number of active (pseudo-)spin degrees of 
freedom. 
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